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■^ . Abstract The linearised general conformal field equations in their first and second 

t^J- | order form are used to study the behaviour of the spin-2 zero-rest-mass equation on 

CNj . Minkowski background in the vicinity of space-like infinity. 



1 Introduction 



In [5], a completely novel finite representation of space-like infinity was proposed. 
In this setting, space-like infinity i° "blows up" to a cylinder / and, consequently, 
asymptotically flat initial data become regular near i° and can be prescribed on 
generic space-like Cauchy surfaces, i.e., the use of hyperboloidal hypersurfaces may 
2_? • ultimately be unnecessary [8]. It turns out that, in this representation, Friedrich's 

/q ' general conformal field equations (GCFE) acquire a very simple form and the cylin- 

der becomes a total characteristic of the system in the sense that there are no radial 
derivatives in the equations restricted to /. In addition, the coordinate location of 
^s , null infinity is now known beforehand and it does not have to be determined during 

the evolution. All these features, namely the regularity of the initial data, the a-priori 
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fixed finite coordinate location of null infinity, the fact that in the Minkowski case 
the entire physical space-time can be covered by one computational domain, and 
the extremely simple form of the evolution equations (especially on the cylinder), 
make the general conformal field equations suitable for numerical manipulations. 
However, as expected, the intrinsic system of evolution equations on the cylinder 
degenerates at the interface of the cylinder / with null infinity. In general, the solu- 
tions generate logarithmic singularities at these regions which are expected to travel 
along null infinity and spoil its smoothness, making it impossible to read-off the 
gravitational radiation at J^ + . The way out of this problem is to prescribe initial 
data that respect the regularity conditions proposed in [5] . 

In this short contribution, we use Friedrich's general conformal field equations 
to evolve generic asymptotically flat initial data near space-like infinity. We begin 
our endeavour from the simplest possible case: linearised gravitational fields on a 
Minkowski background. Although simple, this "toy model" encapsulates most of the 
crucial characteristics of the full non-linear system described above. The initial data 
are evolved as close as possible to the ill-behaved regions 7 ± and study the behaviour 
of the numerical solutions there. This procedure is carried out twice by using the 
linearised general conformal field equations in their first and second order form. 
Analytically the two approaches are equivalent, but their numerical implementation 
could very well differ. The latter statement is partly based on the claim made in [6] 
that writing the equations of general relativity as a system of second order PDE's is 
more advantageous numerically than writing them as a system of first order PDE's. 
According to [6], numerical simulations based on second order PDE's have better 
numerical accuracy and avoid the appearance of spurious waves travelling against 
the characteristic curves. 



2 The spin-2 equations 

We use the spin-2 zero-rest-mass equation for a totally symmetric spinor field <Pabcd 
in the 2-spinor formulation, i.e. V A / A <p ABCD = 0, to model linearised gravitational 
fields on a Minkowski background. Taking the components of the above expression, 
decomposing the five independent components <j> k of the spin-2 into harmonic modes 
(l,m), introducing coordinates (t,r,8,(p) and an adapted spin-frame on the cylinder, 
the above spin-2 equation splits into the eight coupled equations (for details see [1]) 

(l-tK')d t <p k +Kd r <p k -(3K'-(5-k)n)<p k = nc k <p k _ u fc=l:4, 
(l+nOAfc->ri,k + (3K / + (*+l)/i)fc = -/ic i fc + i, £ = 0:3, 



where K(r) = rfi(r), ju(0) = l,c k = y/ 1(1 + 1) - (2-k)(l -k), and ' denotes differ- 
entiation with respect to r. Equations (1), when appropriately combined, split into 
three constraint equations and a symmetric hyperbolic (in the domain \t\ < K ) 
system of five evolution equations, see [1]. 
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By differentiating the above spin-2 equation, one can derive (see [4]) the second 
order spin-2 wave equation □ <J)abcd = 0. Decomposing again the components of the 
spin-2 field into harmonic modes and introducing coordinates (t,r,9,(p), the spin-2 
wave equation splits into a hyperbolic (in the domain \t < K" -1 ) system of five wave 
equations 



(1 -rr)4^- K d rr ((> k + 2tKK'd tr ((> k + 2rKn'd r ((> k 

+ 2[(2 -k)K'-t(K a + r $ 'k' - l -KK")]d t <t> k 

+ [(2-k){KK"+(l-k)K ,2 )+k(5-k)r 2 fl' 2 }fa 

= -LL 2 c 2 k <j> k -rHn'((4--k)c k <j> k+1 +kc k _ l <j> k _ l ). 



(2) 



Notice that the domain of hyperbolicity for both sets of equations is the same, 
namely \t\ < K^ 1 . Thus, as expected, in both approaches the equations degenerate at 
(t,r) = (±1,0), i.e. at the regions / where null like infinity J? ± meets the cylinder. 
In order to compute the same solution as with the first order system (1) we have 
to use all the available information from the first-order system to determine initial 
and boundary conditions for the second order system (2). In [4] it was shown that 
the two formulations are equivalent provided the initial and boundary data for the 
second order system are determined from the first order system. 



3 Numerical results 

The PDE systems (1), (2) are discretised according to the method of lines. An 
equidistant grid on the computational domain D = [0, 1] is used to discretise the spa- 
tial coordinate r. The spatial derivatives are approximated by summation by parts 
(SBP) finite difference operators [7]. The boundary conditions are implemented with 
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1st order 


2nd order 


1st order 


2nd order 


50 


-25.2218 


-27.6006 


-11.1643 


-12.3418 


100 


-29.3956 


-31.5941 


-13.9924 


-15.1743 


200 


-33.7109 


-35.6075 


-16.9978 


-18.1782 


400 


-38.1000 


-39.6068 


-20.1075 


-21.2850 



Table 1 The logarithm of the normalized I 2 norm of the absolute error E, log 2 (| \E\ I2), between 
the exact solution and the solutions computed from the lst-order system (1) and the 2nd-order 
system (2) at time t = I. 



a very simple, but highly efficient, simultaneous approximation term (SAT) penalty 
method [3]. The temporal integration is based on a standard explicit fourth order 
Runge-Kutta scheme. The code has been written form scratch in Python. 



4 Florian Beyer, George Doulis, Jorg Frauendiener, and Ben Whale 

The comparison of the numerical properties of the two approaches is based on 
their ability to reproduce a specific family of exact solutions of the spin-2 equation 
developed in [1]. As was shown in [1,4], the critical sets /+, located at t = 1, can 
be reached in both approaches without loss of the expected 4th order accuracy; in 
addition, the constraint quantities are preserved to sufficient accuracy during the 
evolution. A comparison of the accuracy with which the two approaches numerically 
reproduce the exact solution is shown in Tab. 1 . Better accuracy was achieved in the 
second order formulation, a result that confirms the first claim in [6], The second 
claim made therein that the spurious waves disappear in the second order case is 
also confirmed as the high frequency features disappear in the convergence plots of 
the second order formulation, see [4] . 



4 Conclusion 

In this work, two distinct approaches to the linearised general conformal field 
equations were developed and subsequently implemented numerically. In both ap- 
proaches we managed to reach without loss of accuracy the ill-behaved region /+. It 
is principally not possible to go beyond / ± since the equations loose hyperbolicity. 
A possible way to resolve this problem is presented in [2]. We have also shown that 
the second order formulation of the spin-2 equation leads to a better accuracy by a 
factor of 3-4 and the spurious waves travelling against the characteristics disappear, 
confirming the claims made in [6]. 
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